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Non-SI to SI Units of 
Measurement 


Non-SI units of measurement used in this report can be converted to SI units 
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cubic yards 
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Introduction 


The Earthquake Engineering Research Program embodies several work units, 
one of which is entitled “Geophysical Methods for Site Characterization and 
Measurement of In Situ Properties.” Its objective is to develop geophysical pro¬ 
cedures that provide broader, more economic site coverage while determining 
more engineering parameters than before possible. In previous years, work has 
concentrated on developing techniques for use downstream of dams. However, 
safe, cost-effective design of U.S. Army Corps facilities has been limited by lack 
of accurate information on upstream site conditions. The need exists for devel¬ 
opment of in situ techniques for more accurately determining upstream site con¬ 
ditions relating to earthquake response. In view of the fact that downstream data 
are incorrectly extrapolated upstream, new innovative waterborne geophysical 
techniques must be developed to alleviate this problem. An acoustic impedance 
technique reported by McGee, Ballard, and Caulfield (1995), developed under 
the Dredging Research Program, appears ripe for further development for this 
application. A technically significant upgrade to this technique is the implemen¬ 
tation of a Biot theory approach in assessing acoustic-sediment interactions in 
porous, saturated media. Ramifications of the Biot theory may aid engineers in 
accurately assessing liquefaction potential of Corps structures. This pilot study 
program has been conducted to evaluate new techniques to quantitatively assess 
sediment conditions upstream of man-made dams. 


Site Selection 

Arkabutla Dam and reservoir was selected as the site to demonstrate the fea¬ 
sibility of upstream site characterization using waterborne acoustic techniques for 
a variety of reasons: proximity to the New Madrid seismic zone, existing site 
documentation (borings in the reservoir), support from U.S. Army Engineer Dis¬ 
trict, Vicksburg, and mobilization cost. Further, Arkabutla Dam is the subject of 
a seismic stability analysis. 


General History of Arkabutla Lake 


Arkabutla Lake Project is located in northwestern Mississippi on the Cold- 
water River, a tributary of the Yazoo River. The dam is one of four flood control 
dams in the Yazoo River Basin, Mississippi. The dam is of earthfill construction. 
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approximately 10,700 ft long with a crown width of 40 ft. Approximately 
4,460,000 cu yd of embankment materials were used in constructing the dam. 
Slopes are protected by riprap. Maximum height above the valley floor is 67 ft. 
The outlet works is located near the south abutment of the main dam. It consists 
of a concrete approach, three-gated control structure, transition, single-barrel 
conduit, chute, and stilling basin. Further detailed information, if desired, is 
available from the Vicksburg District. 
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2 Approach 


The methodology used to quantitatively assess the characteristics of the sedi¬ 
ments within Lake Arkabutla is a modified seismic reflection technique that 
relates the engineering properties of sediments to acoustic impedance by pre¬ 
cisely determining the reflection coefficient at each reflection horizon. To 
accomplish this task, a very rigorous Biot-based geoacoustic modeling program 
allowing for detailed investigation of up to thirteen physical fluid/sediment prop¬ 
erties is applied to the acoustic data. The Biot theory and how it is specifically 
applied within the Lake Arkabutla program is provided herein. 


Acoustic Subbottom Profile Data 

Acoustic subbottom profile data (SBP) were collected at 1- and 4-kHz fre¬ 
quencies to map sediment conditions throughout the survey. Figure 1 is an 
example of subbottom imagery from Lake Arkabutla. It graphically shows the 
capability of subbottom profiling in providing stratigraphic and lithological 
assessment of bottom and subbottom sediment environments. 

As the transmitted acoustic signal travels downward through the water col¬ 
umn, reflections are generated at boundaries between two layers of differing 
material properties. In Figure 1, reflections are detected at the lakebed and in the 
subbottom at stratigraphic interfaces. This reflection record is actually from the 
approach channel to the dam outlet works. It shows recent sediment accumula¬ 
tion since lake formation as well as relic streambed sedimentation represented by 
the deeper reflection horizons. 

A more typical subbottom image from Lake Arkabutla, Figure 2, shows an 
east-west survey line starting above the riprap blanket at the toe of the dam. Sur¬ 
veying away from the dam, a layer of unconsolidated sediment overlays an 
acoustically impenetrable horizon probably comprised of the native soil prior to 
dam construction. 

These data are geocoded with accurate position data allowing for delineation 
of the horizontal and vertical extents of all unique sediment units. Another fea¬ 
ture of these data is the contrasting colors (yellow to red color ramp) representing 
relative echo intensity. The higher the contrast, the greater the change in material 
properties. For example, a hard bottom would cause a high intensity (red to 
black) echo whereas the echo from a soft muddy bottom would be a low intensity 
(yellow hue). After calibration of the acoustic systems, these same data are 
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processed analytically to derive quantitative sediment properties such as density 
or porosity. 


Reflectivity and Acoustic Impedance 

Acoustic impedance, Z, is defined as the product of the seismic transmission 
velocity and the density of a material and basically represents the influence of a 
medium’s characteristics on reflected and transmitted waves. Many geotechnical 
properties such as porosity, density, mean grain size, bulk modulus, etc., exhibit 
excellent correlation with impedance. It is possible, therefore, with data of a suf¬ 
ficient signal-to-noise (S/N) ratio, to predict geotechnical properties from normal 
reflectivity data through calculations of the sediment acoustic impedance. For 
example, acoustic impedance of the bottom surface sediment unit is related to the 
reflection coefficient derived at the water-sediment interface through the 
Zoeppritz equation (Zoeppritz 1919) 


Z 2 — z, 


( 1 + R) 

(1 -R) 


( 1 ) 


where 

Z 2 = impedance of second layer (the bottom sediment surface) 

Z; = impedance of the overlying layer (in this case the water column) 
R = reflection coefficient at water/sediment interface. 


Since the sound velocity and density of water are known, Z; can be readily 
calculated. Therefore, to assess bottom impedance (Z 2 ) one must first accurately 
determine the reflection coefficient ( R ) at the water-sediment interface. This is 
computed using what is referred to in underwater acoustics as the sonar equation. 
The sonar equation, discussed thoroughly by Urick (1983), describes the quanti¬ 
tative effects on sonar equipment created by many phenomena peculiar to under¬ 
water sound production. This equation is a design and prediction tool for 
underwater sound applications and relates the effects of the medium, target, and 
equipment. The general sonar equation is given as follows: 

s r =sl-n w - n„+DI+BL+N a (2) 

where 

S R = bottom reflection energy at receiver, db 
SL = total energy of source, db 

N w = 20 x logio (range, meters), db (transmission loss due to spherical 
spreading along the path of propagation) 

Nhyd = receiver sensitivity, db 

Na = amplifier gain, db 
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DI = directivity index of receiving array, db (function of transducer 
beam pattern) 

BL = bottom loss, db = 20 logio(R) 

R = reflection coefficient 

Bottom loss, BL, is evaluated by rearranging Equation 2: 

BL = S r + -SL + N w -N a -DI (3) 

Since all terms on the right side of the equation are either measured directly 
or determined through precise equipment calibrations, BL, and therefore the sur¬ 
face reflection coefficient (BL = 20 log 10 R ) and acoustic impedance (Equa¬ 
tion 1), can be determined. If the desired result is an assessment of the bottom 
surface characteristics, the acoustic solution is complete. All that remains is the 
correlation between the acoustic parameters and physical sediment properties 
through geoacoustic relationships. This correlation is accomplished using a Biot- 
based geoacoustic modeling program. 

For Lake Arkabutla, BL is calculated on all 4-kHz reflection data providing 
continuous survey line coverage of lakebed sediment conditions. 


Biot-based Geoacoustic Modeling (BBSS) 

The Biot theory (1956) was developed to explain acoustic behavior of satu¬ 
rated sedimentary materials accumulating on the seafloor. The applicability of 
the Biot theory to surficial marine sediments has been demonstrated by Stoll and 
Bryan (1970). Ogushwitz (1985) further assessed the applicability of the theory 
to the entire suite of sedimentary materials, from surficial materials with porosi¬ 
ties as high as 90 percent, to well-consolidated materials with porosities as low as 
1 percent. 

Saturated sediment consists of a porous assemblage of sedi m ent grains (the 
“skeletal frame”), whose interconnected pores are filled with water or gas (the 
“pore fluid”). Biot (1956) devised a theoretical model to describe the acoustic 
behavior in such a material. The Biot model treats both the individual and 
coupled behavior of the frame and pore fluid. Energy loss is considered to be 
caused by the inelasticity of the skeletal frame and by the viscosity of the pore 
fluid as it moves relative to the frame. The model predicts that sound velocity 
and attenuation in sediment will depend on frequency, the elastic properties of 
the sediment grains and pore fluid, material porosity, mean grain size, permeabil¬ 
ity, and effective stress. 

Evans-Hamilton, Inc. (EHI) and Ogushwitz developed the Biot-Based Syn¬ 
thetic Seismogram (BBSS) program for conducting Biot-based geoacoustic mod¬ 
eling. The purpose of this program is to correlate measured acoustic responses 
with corresponding sediment properties. The BBSS program specifies values for 
13 physical parameters describing the fluid, grains, and frame properties listed in 
Table 1. Ogushwitz’s (1985) paper provides detailed guidelines for specifying 
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Table 1 

Biot Model Input Parameters 

Fluid Properties 

Density of pore fluid 

g/cm 3 

Bulk modulus of pore fluid 

dyne/cm 2 

Viscosity of pore fluid 

P 

| Grain Properties II 

Density of sediment grains 

g/cm 3 

Bulk modulus of grains 

dyne/cm 2 

Shear modulus of grains 

dyne/cm 2 

Mean grain size 

Cm 

I Frame Properties H 

Porosity 

Dimensionless | 

Dynamic permeability 

cm/sec 

Structure factor 

Dimensionless 

Pore size parameter 

Cm 

Complex bulk modulus of frame 

dyne/cm 2 

Complex shear modulus of frame 

dyne/cm 2 


values of these parameters within the BBSS program. EHI applied the BBSS 
program to the 4-khz acoustic reflection data from Lake Arkabutla, in conjunc¬ 
tion with inferred sediment property information, to develop the geoacoustic 
parameters necessary to create synthetic seismograms based on the specific 
outputs of the BBSS program listed in Table 2. A conceptual schematic of this 
process is presented as Figure 3. 


Table 2 

Biot Model Acoustic Property Outputs 

Properties 

Units 

Wave speed 

M/SEC 

Absorption 

DB/M 

Impedance 

MKS 

Reflection Coefficient 

None 

Signal Amplitude 

Volts 

Phase 

Degrees 


Assessment of Sediment Properties within BBSS 

As shown in Figure 3, acoustic impedance, Z, has been related to sediment 
bulk density. Within the skeletal frame parameters, porosity and the frame 
modulii (Table 1) are the dominant properties affecting wave propagation in 
water-saturated sediments (Ogushwitz 1985). Estimates of the fluid and grain 
properties were arrived at based on standard water properties and from limited 
core data, respectively. The elastic constants and porosity values were then 
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iterated to produce synthetic seismograms matching the actual recorded data. A 
relationship was then developed between measured acoustic impedance and bulk 
density based on the relation between porosity and specific gravity. The density 
function shown in Figure 3 is the same relationship used for the Lake Arkabutla 
study. 
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3 Survey and Equipment 


Survey 

A 1000-ft swath of high resolution subbottom profile (SBP) data, sidescan 
sonar (SSS) data, and high frequency bathymetric data were collected along the 
entire north/south (N/S) length of the upstream side of Arkabutla Dam. Figures 4 
and 5 present actual survey track lines for the 1- and 4-kHz geophysical systems, 
respectively. 4-kHz SBP data were collected along three parallel N/S transects 
approximately 500 ft apart, and thirty-three East/West transects at approximately 
200 ft apart. This provides a balanced matrix of data for producing an accurate 
spatial representation of lakebed sediments. Additionally, 1-kHz SBP data were 
collected along the two easternmost survey lines. SSS data were collected along 
the three N/S 4-kHz tracklines to continuously map lakebed features and struc¬ 
tural appurtenances to the dam. 

The survey was conducted on March 29 through 31. Weather conditions 
were rain with brisk wind producing 0- to 1-ft wind chop on the reservoir. 


Instrumentation 

Due to acoustic interference between systems, the subbottom and SSS sur¬ 
veys were conducted independently. All geophysical data were directly inter¬ 
faced with the survey vessel’s positioning system. Specific instrumentation is 
listed in Table 3. A brief description of each instrument is provided below. 


Survey vessel 

The survey was conducted aboard the survey vessel Nancy Kay owned and 
operated by DIMCO, Inc., Vicksburg, MS. The survey vessel is shown as 
mobilized for the Arkabutla Lake survey in Figure 6. 


Navigation 

A Starlink Differential Global Positioning System (DGPS) was used for 
positioning. Real-time differential corrections were obtained from nearest Coast 
Guard beacon installation, improving horizontal positional accuracy to less than 
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Table 3 

Instrumentation 

System 

Model 

Comments 

Navigation 

Starfink GPS 

Differential GPS with approximately 
±2 m accuracy 

Navigation Software 

HYPACK 

Position logging and survey control for 
all systems 

High-frequency 

echosounder 

Ross Laboratories 
Smartsounder 

200 kHz 

Sidescan sonar 

EG&G SMS-260TH 
W/272-TD Tow fish 

100- and 500-kHz dual frequency sonar 
system 

High-resolution 
subbottom profiling 
(SBP) system 

Datasonics SBP-5000 

4- and 7-kHz operations 

Digital Acquisition 

System 

Triton Isis Data Acquisition 
and Processing System 

Dual processor high-speed acquisition 
system used to collect all geophysical 
data. 


±2 m RMS. Survey navigation, control, and data acquisition was accomplished 
using the HYPACK surveying package by Coastal Oceanographies. This system 
received data from the DGPS and fathometer, performed the appropriate geodetic 
transformations, and then transmitted corrected position and depth information to 
other instrument packages, such as the Isis system, aboard the Nancy Kay. Coor¬ 
dinates used for this investigation are UTM, Zone 15, and are in units of survey 
feet. 


Subbottom profiler 

The higher-frequency acoustic subbottom reflection records were generated 
using a Datasonics SBP-5000 high-resolution “pinger” system. The system was 
operated at 4-kHz to optimize for resolution and depth penetration beneath the 
upper lithological unit. To identify the acoustic data basement, a lower fre¬ 
quency acoustic system referred to as a ‘boomer’ operated below 1-kHz was used 
concurrently. These low-frequency data have poor vertical resolution but allow 
for greater depth of penetration than the pinger system. No statistical data analy- 
is was performed on the boomer data, however, the data proved critical in verify- 
ng data basement. 


Sidescan sonar 

An Edgetech Model SMS-260TH image correcting sidescan sonar system 
with a dual frequency (100- and 500-kHz) transducer tow fish was used for all 
sidescan operations. The system was operated at 500-kHz with a range of 75 m 
per side. Figure 7 shows the Model 272-TD dual frequency tow fish being 
deployed over the bow of the Nancy Kay. Sonar data were provided directly to 
the shipboard Edgetech graphic recorder and in analog form to the ISIS digital 
acquisition system for real-time quality control, display, and data storage. Fig- 
re 8 shows all acquisition system components aboard the Nancy Kay. 
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Digital data acquisition system 

The ISIS shipboard data acquisition and image processing system was used 
to acquire, store, and process all subbottom and side scan data. The Isis system 
also acts as an interface between the navigation and geophysical systems provid- 
ng real-time geocoding of all data. Multiple channels of data were recorded with 
16-bit analog to digital converters operating through multiple sessions of Isis. 
The system was capable of acquisition rates of nearly 33 kHz per channel. 
Graphic displays (see Figure 8) provided real-time quality control of all data 
during acquisition. 
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4 Data Processing and 
Mapping 


Acoustic Reflection Data Records 

Appendix A contains a copy records of the 1- and 4-kHz subbottom profile 
data. The survey track lines in Figures 4 and 5 include the filenames and selected 
ping numbers to correlate geographic location with the data plots. Filenames are 
shown where they begin (at ping 0) with ping numbers annotated every 
400 pings. Shown along with the filename is the CDROM archive disk number 
on which the actual data resides (i.e., D01). Survey direction is in order of 
increasing ping number. 

As stated in Chapter 3, subbottom data were collected at 4-kHz for the entire 
survey for the purpose of characterizing the uppermost sediment. To verify data 
basement, a ‘boomer’ system was operate below 1-kHz along the easternmost 
north-south survey line (Figure 4). For most of the survey, acoustic penetration 
was limited to the bottom of the uppermost sediment unit. The confining layer is 
inferred to be the native soil that was exposed prior to flooding of the reservoir. 
The boomer data did not achieve any more penetration than the 4-kHz system. 
Figure 9 shows a comparison between these two data sets. At both frequencies 
the same reflection horizon confines acoustic penetration. The only location in 
the survey area where considerable acoustic penetration was obtained (Figure 1) 
was within the banks of the approach channel to the intake structure. 


Bathymetry 

Data from the 200-kHz fathometer was used to create a bathymetric surface 
grid and contour map of the survey area. The bathymetry map (Figure 10) shows 
a relatively flat surface. The most significant relief variation is between 
14100500N and 14101000N near the intake tower. The 20- by 20-ft resolution 
digital bathymetric grid was used to create isometric surfaces, sediment thickness 
maps and for the overlays of sediment property data and the surface feature 
maps. 
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Side Scan Sonar Mosaic 


The digital SSS data were processed to create a 0.2-m resolution geocor- 
rected digital mosaic of the entire survey area. A small version of the mosaic is 
presented here as Figure 11. The SSS data supports the subbottom findings 
showing a remarkably uniform, low backscatter reflection surface indicative of a 
soft, smooth, featureless lakebed. At places the more competent native materials 
extrude through the recent sediment deposits as shown by darker reflection 
patches on Figure 11. The mosaic was useful in mapping the extents of the sub¬ 
merged riprap blanket covering the upstream face of the dam below the water 
surface. Any future surveys should include full bottom SSS imagery to assess 
the spatial distribution of surficial sediment conditions. 

Figure 12 shows a close-up view of the intake appurtenances, i.e., wing 
walls, intake passages, and the bottom features in the near vicinity of the 
structure. 

The mosaic is available as a GeoTiff image importable directly into Arclnfo 
GIS. No detailed interpretation of the side scan imagery was performed as part 
of this investigation. 
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5 Geoacoustic Modeling 


Using calibration procedures for data with high S/N ratios, the 4-kHz acous¬ 
tic reflection data were processed to provide estimates of the porosity, density, 
and other associated properties of the sediments found in Lake Arkabutla. Cali¬ 
brations were performed to precisely assess operational parameters of the acous¬ 
tic equipment. Then, in conjunction with known sediment properties, the 
geoacoustic parameters needed to relate acoustic impedance with porosity and 
density were developed. The calibrations procedures are described briefly below. 


Equipment Calibration: Sources and Receivers 

Equipment calibrations involve direct measurements to determine the precise 
transmit and receiving transducer sensitivities. Using the sonar equation (Equa- 
ion 2), determination of transmitter response and total energy at the source, 
source level (SL ), and receiving array sensitivity, , was accomplished by 
rearranging the equation to solve for SL. The specific field procedures and data 
analysis methodology followed precisely as described in McGee et al. (1995). 
Table 4 lists specific equipment parameters measured for the 4-kHz acoustic 
reflection system. These values were used as input parameters to the BBSS geo- 
coustic-modeling program for scaling of output synthetic seismograms and to the 
bottom sediment analysis program (BSED51) used for inversion of the acoustic 
reflection data to sediment properties. 


Table 4 

Equipment Parameters: SBP-5000 

Frequency, kHz 

4 

Output Attenuation, dB 

-10 

SIU Attenuation Setting 

2 ! 

Pulse Length, ms 

0.2 

Wavelet Sample Bin Size 

22 

Source Level, SL, dB re 1 dyne/cm 2 

100 

Receiver Sensitivity, Nhvdr, dB re 1 dyne/cm 2 

-72 
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Determination of Bottom Loss and Surface 
Reflection Coefficient 


Characterization of the sediment surface begins by evaluating the acoustic 
bottom loss, BL, according to Equation 3. By substituting into Equation 3 the 
equipment parameters from Table 4 (SL and N^, N w , and DI), and evaluating 
the surface data throughout the survey area, a contiguous assessment of sediment 
acoustic response is accomplished. The acoustic solution for the bottom sedi¬ 
ments is now complete. For the surface data, the remaining task involves cor¬ 
relation of acoustic response to sediment properties, hence establishing the 
geoacoustic calibration. All geoacoustic assessments are primarily related to the 
acoustic impedance, Z 2 . 


Biot Modeling 

Biot parameters 

As stated in Chapter 2, EHI developed the BBSS program for conducting 
Biot-based geoacoustic modeling to correlate measured acoustic responses with 
in situ sediment properties. These properties are typically confirmed through 
sediment core analysis. Specific implementation of the program is through a 
framework of forward modeling requiring correlation of predicted acoustic 
responses based on specified sediment conditions with actual measured acoustic 
returns. Known or inferred sediment parameters (Table 1) are used as inputs to 
the model and the desired quantities are manipulated to give a matching seismic 
response to the actual acoustic data (Table 2). Verification is provided by coring 
and sample analysis. For the Arkabutla study, very limited sediment data were 
available. Core data, included in Appendix B, provided soil and sediment prop¬ 
erties in terms of the Unified Soil Classification System, i.e., soil type, Atterburg 
limits, moisture contents, etc., but no laboratory measurements of any of the 
thirteen physical Biot parameters. Although quantitative property information 
was not available, the information did provide general constraints within which 
the Biot theory could be used in a predictive mode. 


Sediment categorization 

Table 5 presents the basic categories used to report sediment conditions in 
Lake Arkabutla. Sediment categories are grouped according to a density range 
associated with sediment types. These definitions were developed by correlating 
acoustic response characteristics with the available core data through the BBSS 
program. Given the lack of sample data, sediment bulk density is considered the 
most reliable sediment parameter to predict. Sediment density is closely related 
to the elasticity of the sediments and can be calculated from acoustic impedance; 
a parameter that represents the influence of the medium’s characteristics on 
reflected and transmitted waves. For all categories density is the most accurate to 
derive acoustically and best represents insitu conditions. 
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Table 5 

Sediment Descript 

ion 

Density, g/cm 3 

Basic Sediment Description j 

<1.2 

Fluidized mud, very soft 

1.2 -1.4 

Unconsolidated mud, saturated clay, silt 

1.4-1.7 

Clay-silt, clay-silt-sand mixtures 

1.7-2.0 

Stiff clayey sand, silty sand, fine sand 

2.0 - 2.2 

Dense clays, unsaturated sediments, possible organics 

N/A 

Riprap blanket 


For this discussion, fluid mud is defined as sediments with density values 
less than 1.2 g/cm 3 . Sediments below 1.2 g/cm 3 possess physical properties typi¬ 
cal of fluidized unconsolidated sediments. These sediments were modeled with 
porosities greater than 0.85 and negligible frame and shear modulii (specific 
modeling results are discussed later. Sediments in the 1.2 -1.4 g/cm 3 range 
responded geoacoustically as an inelastic frame (reasonable frame modulii) and 
have physical properties representative of freshwater clays and silts. 

The sediments with densities higher than 2 g/cm 3 are potentially dense clays 
or even unsaturated sediments present before the filling of the reservoir. As will 
be shown in the next section, acoustic responses from these sediments were often 
outside the response range typical of marine sediments. 


Modeling results 

After reviewing all the data, three sites were chosen to represent the different 
acoustic environments found underneath Lake Arkabutla. These sites were 
selected to demonstrate the ability of the technique to meet the objectives of this 
investigation as well as its potential limitations. 

Site Arkl6a. This site is located immediately upstream (east) of the intake 
tower within the confines of a submerged entrance channel along survey file 
AR040016 in Figure 5 (refer also to Figure 24 for actual location). Figure 1 is 
the actual 4-kHz subbottom profile data from which the model site was selected. 
The model site was chosen near the area of deepest acoustic penetration. As can 
be seen by this reflection profile (Figures 1 and 24), nearly 17 ft of acoustic pene¬ 
tration was achieved at 4-kHz. It is likely this feature was part of the relic chan¬ 
nel where erosion and depositional cycles occurred prior to filling the reservoir. 
An actual acoustic signal representative of this depositional environment is 
shown in Figure 13. 

A five-layer model (water column plus four sediment layers) was required to 
match the data. An acceptable calibration (Figure 13), requires that actual and 
predicted signal amplitude and phase match accordingly. The synthetic seismo¬ 
gram (red line in Figure 13) is produced after convolution of a field-measured 
source wavelet (accomplished during equipment calibration) with the reflection 
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sequence modeled according to Biot theory. Table 6 presents the specific input 
parameters required for the BBSS model. 


Table 6 

Biot Model Geoacoustic Input Parameters; Site Ark16a 8 


niviAM 

|* Ml 


■H'.'-IIEB 


Label 

Water 

Soft Mud 

Clay-silt- 

sand 

EES® 

Consolidated | 
Sediment | 

Travel Time to bottom of layer, 
ms 

12.20 

13.30 

16.15 

18.54 

N/A 

1 Fluid 1 

Fluid Density, g/cm 3 

1.0 

1.0 

1.0 

1.0 

1,0 1 

Fluid Bulk Modulus, dynes/cm 2 , 
xE+11 

0.215 

0.215 

0.215 

0.215 

0.215 | 

Fluid Viscosity 

0.019 

0.019 

0.019 

0.019 

0.019 I 

i Grain | 



2.53 

2.65 

2.65 

2.65 j 

Grain Bulk Modulus (Real), 
xE+11 


5.40 

4.50 

4.50 

3.64 | 

Grain Shear Modulus (Real), 
xE+11 


2.70 

2.10 

2.10 

1.80 

Mean Grain Size, cm 


0.0016 

0.0062 

0.016 

0.016 J 


Frame 




Porosity 


0.84 

0.55 

0.47 

0.36 j 

Permeability, cm/sec, xE-05 


0.04 

0.003 

0.0007 

0.00006 I 



No 

Frame 

2.50 

2.50 

102.0 

Frame Bulk Modulus (Imag), 
dynes/cm 2 , xE+09 


No 

Frame 

0.027 

0.027 

1.12 

Frame Shear Modulus (Real), 
dynes/cm 2 , xE+09 



1.15 

1.15 

47.0 

Frame Shear Modulus (Imag), 
dynes/cm 2 , xE+09 


No 

Frame 

0.037 

0.037 

1.50 


Sediment layers range from soft muddy sediments to a dense, high imped¬ 
ance layer considered the acoustic ‘basement.’ The soft mud blankets the lake- 
bed and is found throughout the survey area. This mud represents recent 
deposition and shows properties typical to unconsolidated fine-grained sedi¬ 
ments. This layer was modeled with negligible frame and shear bulk modulli 
which account for the lack of shear strength normally found in consolidated sed¬ 
iment materials. Porosity was determined to be 0.84 and wave velocity and 
acoustic impedance 1429 m/sec (indicates near saturation) and 1778 mks, 
respectively. 

Underlying this soft mud layer is a unit of clay/silt/sand (Table 6, Layer 3) 
with a porosity of 0.55. This results in a high-impedance contrast with the upper 
layer. Wave velocities increased to 1539 m/sec. Average bulk density modeled 
at 1.74 g/cm 3 . 

Beneath the clay/silt/sand unit is a more competent layer (Table 6, Layer 4). 
An increase in density from 1.74 g/cm 3 to 1.87 g/cm 3 could result from a higher 
percentage of sand in the material. 
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A high impedance layer (> 6300 inks), represented by the high-amplitude 
reflector in Figure 13, is the confining acoustic interface and was modeled as a 
consolidated sediment with very high frame bulk and shear modulli (Layer 5 in 
Table 6). Polarity of the reflection coefficient remained positive indicating an 
increasing impedance environment (i.e., a harder material than the overburden 
sediment). No penetration was achieved below this layer with any of the acoustic 
systems. 

The sediment cores did not provide enough information to establish absolute 
verification of the model. However, based on past experience and general guid¬ 
ance regarding the likely sediment environment at Lake Arkabutla, the results 
should be considered reasonable estimates of actual conditions. Except for 
pockets of trapped gas (i.e., organics) and consolidated, unsaturated sediments, 
estimates of the acoustic responses (impedance, velocity, phase, and reflection 
coefficient) should be accurate. 

Site Arkl4a. This site represents a potential difficulty in using normal- 
incidence reflection techniques where the sediments were previously exposed 
native soils submerged after filling of a reservoir. Figure 2 from Chapter 1 is the 
profile data from line 14 and shows a soft mud layer overlying a highly reflective 
and acoustically impenetrable layer at Model site Arkl4a (actual location shown 
in Figure 23). Figure 14 shows the upper sediment unit modeling as a soft mud 
identical to the uppermost unit model at site. Acoustic penetration is arrested at 
the bottom of this soft mud unit. Specific BBSS input provided in Table 7 does 
not list any sediment data for Layer 3. Attempts to model the actual data failed 
when using saturated sediment parameters. The matching response shown in 
Figure 14 was obtained by establishing the bottom boundary of the model at this 
interface and setting the reflection coefficient to -0.80. The large negative reflec¬ 
tion coefficient indicates a decreasing impedance sequence and a drastic change 
in sediment type. Normally, this would be caused by gas bubbles trapped in the 
sediments. However, the fact that this area was most likely exposed pasture 
lands prior to flooding the reservoir, the acoustic response could indicate very 
low moisture content clays similar to hardpan. Either of these conditions could 
preclude one’s ability to achieve sufficient penetration into these materials, par¬ 
ticularly at 4-kHz. Also, this modeling technique is only valid in gas-free satu¬ 
rated sediments and could be a limiting factor at other sites regardless of seismic 
frequency. 

Site Arkl4b. This site is located approximately 200 ft past the end of the 
riprap blanket covering the toe of the dam (Figure 2 and Figure 23) along a short 
shelf that is possibly an engineered extension of the upstream face of the dam. 
These data indicate a more fully saturated sediment environment conducive to 
modeling hydroacoustic data. There are no measured subbottom reflectors in the 
vicinity of site Arkl4b (Figure 2). Model comparison is presented in Figure 15 
and the input parameters are listed in Table 8. The sediment models as a stiff 
clay/silt material with a porosity of 0.58, density of 1.64 g/cm 3 , and reasonable 
frame bulk and shear modulli. The data are ‘in-phase’ with the source wavelet 
indicating a ‘gas-free’ environment. The core data (Appendix B) were taken in 
the vicinity of this site, i.e., near the toe of the dam, and show sands and clays 
with relatively low blow counts supporting the model predictions. 
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Table 7 

Biot Model Geoacoustic 

Input Parameters; Site Ark14a I 


Layer 1 

Layer 2 

Layer 3 

Layer 4 

Layer 5 

Label 

Water 

Soft 

Mud 

Basement: 

Organic/ 

Unsaturated 

Not 

Used 

Not 

Used 

Travel Time to bottom of layer, ms 

9.13 

9.98 




II Fluid J 

Fluid Density, g/cm 3 

1.0 

1.0 




Fluid Bulk Modulus, dynes/cm^, 
xE+11 

0.215 

0.215 




Fluid Viscosity 

0.019 

0.019 




1 Grain II 

Grain Density, dynes/cm 2 


2.53 




Grain Bulk Modulus (Real), xE+11 


5.40 




Grain Shear Modulus (Real), 
xE+11 


2.70 




Mean Grain Size, cm 


0.0016 




|| Frame H 

Porosity 


0.84 




Permeability, cm/sec, xE-05 


0.04 




Frame Bulk Modulus (Real), 
dynes/cm 2 , xE+09 


No 

Frame 




Frame Bulk Modulus (Imag), 
dynes/cm 2 , xE+09 


No 

Frame 




Frame Shear Modulus (Real), 
dynes/cm 2 , xE+09 


No 

Frame 




Frame Shear Modulus (Imag), 
dynes/cm 2 , xE+09 


No 

Frame 





Table 8 

Blot Model Geoacoustic Input Parameters; Site Ark14b 



Layer 1 

Layer 2 

Layer 3 

Layer 4 

Layer 5 

Label 

Water 

Stiff 

Clay/Silt 

Not 

Used 

Not 

Used 

Not Used 

Travel Time to bottom of layer, ms 

7.47 

N/A 




I Fluid _i 

Fluid Density, g/cm 3 

1.0 

1.0 



| 

Fluid Bulk Modulus, dynes/cm*, 
xE+11 

0.215 

0.215 



I 

Fluid Viscosity 

0.019 

0.019 



u 

1 Grain II 

Grain Density, dynes/cm 2 


2.65 




Grain Bulk Modulus (Real), xE+11 


5.40 




Grain Shear Modulus (Real), 
xE+11 


2.70 




Mean Grain Size, cm 


0.0016 




1 Frame II 

Porosity 


0.58 




Permeability, cm/sec, xE-05 


0.003 




Frame Bulk Modulus (Real), 
dynes/cm 2 , xE+09 


2.50 




Frame Bulk Modulus (Imag), 
dynes/cm 2 , xE+09 


0.027 




Frame Shear Modulus (Real), 
dynes/cm 2 , xE+09 


1.15 




Frame Shear Modulus (Imag), 
dynes/cm 2 , xE+09 


0.037 
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Modeling summary 


Three sites were modeled to describe general sediment conditions within 
Lake Arkabutla. Insufficient core data were available to intensively evaluate the 
BBSS program and to provide verification of predicted sediment types. Acoustic 
penetration was limited with the subbottom profiling systems used. Future sur¬ 
vey programs should include a small sampling/coring program in conjunction 
with the geophysical surveys to calibrate and verify the analysis. 

Given the limited ground truth data, the approach provided a method to rigor¬ 
ously interrogate actual sediment conditions from acoustic data. The Biot theory 
is a proven model. Acoustic calibration methods have been developed and 
employed here to quantify precisely the net acoustic response in a given sediment 
sequence. For data with sufficient signal response (Arkl6a), the method pro¬ 
vided reasonable estimates of sediment characteristics. The method also pro¬ 
vided a critical evaluation of ‘non-standard’ sediment environments, such as site 
Arkl4a, where it was not possible to match actual acoustic data using standard 
sediment parameter inputs. 


Inversion of Acoustic Data 

Acoustic reflection data are processed through seismic inversion techniques 
to produce continuous sets of geopredicted sediment property values. This is 
accomplished through a software package developed by EHI, the BSED51 prog¬ 
ram, designed to incorporate calibrated acoustic data with modeled geoacoustic 
parameters. The BSED51 program was used to process all 4-kHz acoustic data. 

It provides graphics and digital sediment property output files and can be deliv¬ 
ered as input data for all visualization and CADD work. Figure 16 shows density 
output graphic from the BSED51 program for the water-sediment interface along 
survey line 14 (Figure 2). Selected porosity (P) and density (p) values deter¬ 
mined from the Biot calibration have been included for comparison with surface 
densities computed using BSED51. It shows excellent correlation between the 
modeled and processed results. The red dashed-line represents detected depth at 
the reflection interface used in the density calculations. Detection levels are con¬ 
trolled by the ‘Threshold’ variable and is analogous in function to the sensitivity 
control on a standard echosounding system. Processing input parameters are 
listed in the green control panel on the left side of the figure, and includes equip¬ 
ment calibration values (Table 4), detection threshold controls, transmitter/ 
receiver geometry, and statistical processing controls. The geoacoustic calibra¬ 
tion function, listed as ‘Mobile Bay’ in the control panel, is the impedance versus 
density model used for the Arkabutla Lake data. This relationship was derived 
during the geoacoustic-modeling phase of a previous survey performed in Mobile 
Bay, Alabama (McGee 1998). Acoustic response was similar to this survey in 
that a layer of unconsolidated fine sediments, or fluid mud, overlaid more compe¬ 
tent native sediments. All 4-kHz data were processed through the BSED51 prog¬ 
ram to create matrices of position coordinates, depth, acoustic impedance, 
density, and porosity. 
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The BSED51 program was also used to conduct inversion on selected sub¬ 
bottom reflection interfaces. Accounting for two-way transmissivity and absorp¬ 
tion through the upper layer the reflection coefficient and acoustic impedance can 
be calculated for subsequent deeper layers. The geoacoustic calibration function 
is the same as for the water-sediment interface. Figure 17 shows the graphical 
output screen from BSED51 along with selected BBSS reflection coefficient, R, 
value for model site Arkl4a. Density spikes represent acoustic responses outside 
the boundaries established for a natural sediment response. Conditions at these 
locations seem to be unsaturated or consolidated clays and silts. The presence of 
interstitial gasses is also a possibility. 
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Data Presentation 


A brief description of physical conditions is provided to assess the spatial 
distribution of sediments within the surveyed portion of Lake Arkabutla. The 
purpose here is to present a range of possible data products available to assist in 
interpretations of geophysical survey data. A digital mosaic of side scan sonar 
imagery (Figure 11) was created to provide a contiguous picture of surficial sedi¬ 
ment characteristics. Lakebed bathymetry data was compiled into a 20-ft by 
20-ft resolution surface grid for contouring and as a base for overlays of specific 
sediment property information. Output from the BSED51 program was gridded 
for bulk density and sediment thickness to describe the actual sediment condi¬ 
tions. Figure 18 presents selected isometric views of sonar imagery, bathymetry, 
and sediment density with respect to water depth. This figure shows how this 
type information can be combined to provide enhanced interpretations of actual 
bottom conditions. For example, on the right side of Figure 18 is a small hard 
reflector (black contrast) near the front of the image. This corresponds to a small 
spike in the bathymetric data and a small high-density patch in the density image. 


Lakebed Sediment Density and Thickness 

Surficial sediment density was continuously mapped over the extent of the 
survey area (Figures 18 and 19). Density is for the entire thickness of the upper 
sediment unit and shows a uniformly distributed layer of 1.2 -1.4 g/cm 3 material. 
Density increases nearer the toe of the dam (western edge of map). Density is 
not reported for the riprap placed near the toe of the dam. 

For Arkabutla, and other reservoirs, these data offer a complete picture of the 
distribution of sediments accumulating on the bottom of the lake. Combining 
physical property data such as density, with sediment thickness maps, it is possi¬ 
ble to calculate the total volume of mud accumulated. Figure 20 presents a mea¬ 
sured thickness for the upper sediment unit. This map was developed from the 
subbottom profile data by combining reflection interfaces from the digital subsur¬ 
face model. Except for survey line 16 near the intake tower entrance, only a 
single subbottom interface at the bottom of the mud layer was detected. The data 
was not sufficient to develop more elaborate lithological models of Lake 
Arkabutla. 
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Sediment Profiles 


Six representative subbottom profile sections are presented in Figures 21 
through 26. Each figure shows the actual subbottom profile data with the associ¬ 
ated two-dimensional interpretations of lithology. This presentation is the basic 
‘fence diagram’ format familiar to most geophysicists. The sediment units are 
grouped according to density range. After modeling, a basic sediment descrip¬ 
tion was developed for each density range as presented in Table 5 (Chapter 5). 
Had there been confirming core sample data, more precise descriptions would 
have been developed, possibly relating to specific seismic issues. The profiles 
are self-explanatory and will only be briefly discussed following. 

A 2- to 3-ft thick layer of soft mud overlays native soils along the eastern 
two-thirds of the survey area. No penetration was obtained through the native 
soils, resulting in a simple two-layer lithology. The western edge of the profiles 
extends over a riprap blanket along the toe of the dam. This zone of riprap is 
characterized by high-amplitude acoustic returns followed by numerous multiple 
reflections. Acoustic response in native soils tend to be high-amplitude with 
inconsistent wavelet polarity indicative of organic/gassy or unsaturated 
environments. 

As shown by Figures 21 through 26, subbottom profiling allows for continu¬ 
ous surveys of subbottom conditions from which an optimum sediment sampling/ 
coring plan can be developed assuring sampling from all sediment units. 

Through geoacoustic modeling of the sediment and acoustic data, physical prop¬ 
erties can be determined for all survey data providing accurate data over large 
geographic areas. This would reduce or even eliminate the potential of erroneous 
extrapolation of sediment conditions between isolated core locations. 
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7 Conclusions and 
Recom mendations 


Prior to the filling of Lake Arkabutla, the area upstream of the dam was com¬ 
prised of agricultural lands. The only sediments to be found were those depos¬ 
ited within the banks of the streams meandering through the watershed. Except 
for the relic streambed, the former surface soils act as the confining acoustic 
interface (Figure 2), limiting penetration to the bottom of the layer of sediments 
accumulated since filling of the reservoir. Within the meandering streambed, 
nearly 17 ft of acoustic penetration was achieved for the 1- to 4-kHz frequency 
bandwidth chosen. In these areas, sediment characterization was reasonable as 
shown by Figures 1,13 and 24. 

The instruments selected for this initial investigation were not adequate due 
to frequency and source level to overcome attenuation caused by the confining 
layer. This confining layer is likely comprised of compacted, unsaturated soils 
possibly containing decomposing organics found over most of the project site. 
This may be a typical configuration for man-made reservoirs. For an acoustic 
solution, a lower frequency system is recommended. 

Specific conclusions and recommendations for the Lake Arkabutla sediment 
characterization feasibility study are provided following. 


Conclusions 

Biot theory 

The Biot theory is a proven method (see references). For Arkabutla survey 
data with sufficient signal response, the method provided, in the absence of quan¬ 
tified ground truth from the core data, reasonable estimates of sediment density, 
porosity, acoustic wave speed and attenuation were obtained. 

The method of forward modeling to develop site specific transfer functions 
relating acoustic response to sediment properties provided a critical evaluation of 
the sediment environment at Lake Arkabutla. This allowed for analysis of ‘non- 
standard’ sediment configurations. 
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Sediment bulk density, computed from acoustic impedance, is considered the 
most reliable sediment parameter for Arkabutla and best represents in situ 
conditions. 


Core sampling 

No sampling or coring was conducted in conjunction with this investigation. 
Core data retrieved from near the upstream toe of the dam (Appendix B) was 
provided from a previous investigation (reference if needed). 

No laboratory measurement of any of the thirteen physical Biot parameters 
was available from the existing core data. This provided only general guidance 
from which the Biot theory could be used in a predictive mode. 

Core positioning was not coordinated in conjunction with the geophysical 
survey nor was it possible to accurately determine their precise location. 


Acoustic response 

Characterization of sediments accumulated since filling of the reservoir was 
accomplished in terms of density and sediment thickness providing information 
useful in accurately determining the rate and composition of sedimentation in 
Lake Arkabutla. The data were collected rapidly and continuously, providing 
full coverage mapping of sediment conditions within the lake. 

Information regarding the reflection coefficient for the water-sediment inter¬ 
face, attenuation through the upper sediment layer, and reflectivity of the confin¬ 
ing layer was obtained. 

Within the former riverbed, acoustic penetration and signal response was 
sufficient to characterize sediment conditions down to nearly 17ft below the lake 
bottom using the 1- and 4-kHz systems. 

An acoustically confining layer of compacted unsaturated soils limited 
acoustic penetration to the bottom of the sediment layer. This layer did not 
respond acoustically as a standard sediment. Both the 1-kHz and 4-kHz systems 
were unable to penetrate through this native soil unit. 


Recommendations 

A follow-on field study with an acoustic source operating in the seismic 
range, i.e., 100-500-Hz should be conducted to fully evaluate the feasibility of an 
acoustics approach. This would likely involve use of a small airgun source. 

A ground penetrating radar (GPR) system should be included as part of the 
field investigation to verify the presence of acoustic-inhibiting conditions such as 
decomposing organics. 
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A coring program resulting from the geophysical investigation is recom¬ 
mended. Core locations should be selected based on survey results. If good core 
data already exists, the geophysical survey must be designed to collect acoustic 
information precisely over the core location. 

Cores should be analyzed for density, porosity, sound velocity, and specific 
gravity along with a complete grain size analysis. This will provide sufficient 
information for calibration of geoacoustic conditions using the Biot theory. 
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Figure 1. Example subbottom profile from Lake Arkabutla 
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Figure 3. Biot-based modeling concept 















Lake Arkabutla 4-kHz survey track lines. 
















Figure 6. Survey vessel Nancy Kay at Lake Arkabutla 
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Figure 7. Side scan tow fish deployed over bow of Nancy Kay 
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Figure 8. Geophysical data acquisition 





































10. Bathymetry map. Lake Arkabutla. 














;ure 11. Side scan sonar mosaic, Lake Arkabutla. 










Figure 12. Side scan image of approach to intake 






Water 



Consolidated 

sediment 


Vel = 1429 m/s 
lmpedance=1778 mks 
Density=1.24 g/cc 
Porosity=0.84 


Vel = 1539 m/s 
Imped a nce=2682 mks 
“Density=1.74 g/cc 
Porosity=0..55 


Vel = 1583 m/s 
lmpedance=2968 mks 
Density=1.87 g/cc 
Porosity=0.47 



3071 m/s 


lmpedance=6314 mks 
Density=2.06 g/cc 
Porosity=0.36 





























Figure 15. Geoacoustic model: SiteArk14b 
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Figure 16. BSED51 output at water-sediment interface 
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Figure 17. BSED51 output at data basement 
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Selected isometric views. 














19. Surficial sediment density distribution. 
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Sediment profile. Line 21. 







































































Appendix A 

1-kHz Subbottom Profile Data 


Appendix A 1-kHz Subbottom Profile Data 
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